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We present the calculation of the hyperon forward spin polarizability 70 using manifestly Lorentz 
covariant baryon chiral perturbation theory including the intermediate contribution of the spin 3/2 
states. As at the considered order the extraction of 70 is a pure prediction of chiral perturbation 
theory, the obtained values are a good test for this theory. After including explicitly the decuplet 
states, our SU( 2 ) results have a very good agreement with the experimental data and we extend our 
framework to SU( 3 ) to give predictions to the hyperons’ 70 values. Prominent are the E” and 
baryons as their photon transition to the decuplet is forbidden in SU( 3 ) symmetry and therefore 
they are not sensitive to the explicit inclusion of the decuplet in the theory. 
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I Introduction 


From the experimental study of Compton scattering on a baryon target one can extract relevant information about 
the inner structure of baryons. With the help of the sum rules for integral characteristics of the cross sections, very 
important observables like polarizabilities can be assessed. The focus of this work is the forward spin polarizability 
7o, which represents the deformation of a hadron relative to its spin axis when scattering photons in the extreme 
forward direction. It is related to the photo-absorption 7 V cross sections ( 73 / 2 , 1/2 with total helicities 3/2 (for parallel 
photon and target helicities) and 1/2 (for antiparallel photon and target helicities) via the sum rule in Ref. [T] 
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originally found in Ref. [ 2 ]. The energy ujq is the threshold for an associated neutral pion in the intermediate state. 
Experimental results for the proton 70 were obtained in Ref. [3] and, more recently, in Ref. [5]. Furthermore, dispersion 
relation studies have been performed in Ref. [5] for both nucleon spin polarizabilities. 

Concerning the theoretial approach, the nucleon’s structure has been thoroughly studied with the help of effective 
field theories on Compton scattering data in Refs. 0,0 and [8]. The spin-dependent piece of the amplitude 
attracted particular interest. The term proportional to 00 ^ (cu is the photon’s energy) contains the whole information 
about 7o, via the master formula 
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as described in Refs. DP and m- Here a is the vector of Pauli matrices, e and are the polarizations of incoming 
and outgoing photons, respectively, a is the fine-structure constant and e the elementary charge. 

Early calculations in models of chiral perturbation theory (ChPT) that include only nucleonic intermediate states 
have been performed both in a heavy-baryon approach as well as in fully covariant calculations as in Ref. m- In 
Refs. [ 12 ] and [13] the theory was extended such as to include isospin-3/2 intermediate states, namely the A(1232) 
resonance. It was found that the inclusion of the latter state greatly improved the convergence between theory and 
empirical evidence. 

When considering ChPT in SU(3) models, valuable predictions about the hyperons’ polarizabilites can be calculated, 
where there are no experimental data available yet. First results with the help of heavy-baryon ChPT were obtained 
in Ref. [10] and later improved in our work. Ref. [TJ] . The predictions are expected to be more reliable when using a 
fully covariant model and introducing the A(1232). 

Along this line in this work we perform a calculation of the amplitude including corrections induced 

both by a fully covariant version and intermediate spin 3/2 states in a full extension to SU(3) flavor. The leading 
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ChPT order for the quantity 70 is a calculation. When including the A(1232) resonance, the study of Ref. [12] uses 
the so-called small scale expansion introduced in Refs. m and m- We opt for a different power counting scheme, 
following Ref. m and, furthermore, we introduce the couplings in a consistent dynamics, using the full A(1232) 
propagator as in Refs. m, na, m and [ 21 ] . A study comparing different effective field theoretical models has been 
performed in Ref. m 

The outline of this paper is as follows. In Section |TI| we give a theoretical introduction to the appropriate chiral 
Lagrangians and the power-counting used. The kinematical considerations and assumptions for the calculation of 70 
are presented in Section [ITT] and the results are discussed in Section [TV| Finally, we briefly summarize in Section [y] 


II ChPT involving pseudoscalar mesons, baryons and photons 


For the description of hyperon polarizabilities we use a manifestly Lorentz covariant SU(3) version of chiral pertur¬ 
bation theory involving pseudosclar mesons, baryons and photons (see details in Refs. [23] and [23]). The lowest-order 
chiral Lagrangian involving pseudoscalar mesons 0, baryons B and photons reads 

c = cfl + c%, ( 3 ) 

where 

is the O(p^) meson Lagrangian and 

= Tr {B{i 0 - m)B) + ^Tr (^^75 {u^, B}) + ^Tr (^^75 K, B]) ( 5 ) 

is the 0{p^) Lagrangian including baryons. The symbols [ ] and { } occurring in Eq. [^and in the following denote 
the commutator and anticommutator in flavour space, respectively. The vielbein is given by i with 

= U = exp where = d^u — \{y^ + clii)u + — a^) and = d^B + [T^^B] are the covariant 

derivatives acting on meson and baryon octet fields, respectively, and m denotes the baryon octet mass in the chiral 
limit. The chiral connection is given by F^ = , d^u] — + cl^)u — — cl^)u^• Since we are working with 

photon fields, we set to ee^Q and to 0. The constant Fq is the meson decay constant in the chiral limit and the 
low-energy constants D and F are determined from hyperon (3 decays, where the combination F F D corresponds to 
the low-energy constant qa in the SU(2) limit. The explicit form of the 3x3 charge matrix Q, meson (p and baryon 
B matrices is given in Appendix [A] The term Tr(x+) is responsible for the explicit breaking of the chiral symmetry 
due to the finite quark masses 

Tr(x+) = Tr (xt/t + , (6) 

where in our case y = and M is the meson mass. The power-counting scheme followed here gives the order 


N = ANl + J2dNd- 2P^ - Pb 

d=l 


(7) 


to a diagram, where Nl stands for the number of loops, Nd for the number of vertices from Lagrangians of order d, 
and P(j) and Pb for the number of meson and baryon propagators, respectively (see also Ref. [25]). 

In this work we also include isospin-3/2 resonances, which give significant corrections to the full amplitude. The 
relevant terms of the Lagrangians that couple these decuplet fields to the octets of baryons and mesons are partly 
given in Refs. Hg, HZ], HH], where the kinetic term reads 

£7 = F'”' - MaY''] Af 7 (8) 

We added the missing couplings by extending the known vertices from SU( 2 ) (see Refs. [22] and [21]) to SU(3): 
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where are the decuplet states (see details in Appendix]^, is the self-dual stress tensor of the 

electromagnetic field, and the Dirac tensors 7 ^^ and 7 ^^^ are specified in Appendix |a| 

The couplings C and qm are low-energy constants and Ma corresponds to the decuplet mass in the chiral limit. 
Note that the constant C corresponds to the low-energy constant Ha of SU( 2 ) in Refs. [29] and [ 21 ] by the conversion 
C = — . This definition of Ha differs by a factor 2 from the definition found in Ref. m- This low-energy constant 

is extracted from the strong decay of the decuplet into the baryon octet and has been determined to be = 2.85 in 
Ref. [30] for SU( 2 ) and C = —0.85 in Ref. [28] for SU(3). It is also important to mention that the numerical value for 
the coupling constant qm has not yet been studied when extending the model to SU(3). Therefore the quality of the 
predictions very much depends on its correct value. We follow the method of Ref. m and estimate the value of qm 
by calculating the width of the electromagnetic decay of the A(1232): 


t^EM 
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where Im(S™) is the imaginary part of the electromagnetic A(1232) self-energy amplitude. Therefore, using the 
relation r™/(r™ + p^trong^ _ q 55 % q 55 % ^nd the strong decay width ± 2)MeV, we get the value 

Qm = 3.16 ± 0.16. Since data on the electromagnetic decays of the full decuplet are sparse and contain large errors, 
a determination of qm in the SU(3) version is not viable. We therefore also fix qm in SU(3) to the A —^ 7 A decay, 
i.e we also use the value of qm = 3.16 ± 0.16 here. We should keep in mind that the central value of qm will change 
when going from SU( 2 ) to SU(3), but we expect that with the present sizable error on qm this value is included. 

For the covariant derivative we use 


{D^ct>r^ = T _ (^2) 

The factor ^ comes from the definition of the meson-octet matrix. When introducing the decuplet fields into the 
chiral theory, an additional small parameter of the ChPT expansion appears, 6 = Ma —m^ 300MeV. Therefore, the 
counting scheme has to be revised. Here we follow the 5 counting scheme from Ref. HZ], where 5 ‘^ is counted as 0{p). 
It is adequate for the low-energy range close to pion-production threshold. Hence, one obtains the rule 

oo ^ 

N = mL + Y, dNd - 2P^ -Pb- -Pa, (13) 

d=l 


where now Pa is the number of A propagators. Since the A(1232) is a spin-3/2 resonance, it does not have the normal 
Dirac-propagator form for spin-1/2 particles, but takes the Rarita-Schwinger form given by 
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where d is the number of dimensions of the Minkowski space, which after dimensional regularization is set to d = 4. 

The calculations in this work are done up to order p^ for the isospin-1/2 counting scheme, which is the leading-order 
calculation for the 70 observable, and up to order in the isospin-3/2 counting scheme. This choice is due to the 
fact that in the isospin- 1/2 sector the first contributions appear at loop level, which corresponds to the order if 
all the coupling vertices are extracted from the lowest-order Lagrangian. Instead of going to higher-order couplings 
and therefore obtaining loops of order p^, we included the isospin-3/2 sector, which due to its lower order — the 
leading-order diagrams are at order p^/^ — is expected to dominate over those contributions. In addition there is 
the advantage that fewer additional low-energy constants are needed than for the case of higher orders. All the 
constants’ values chosen for this work are given in Table In order to obtain the SU( 2 ) results, we simply put all 
the channels with non-vanishing strangeness to 0 and keep only those channels involving nucleons, pions and the 
isospin-3/2 quadruplet. 


Ill The forward spin polarizability 70 

The diagrams contributing to the forward spin polarizability 70 are shown in Figs. and The calculation of the 
amplitudes corresponding to each of these diagrams is performed in the rest frame of the baryon. To calculate the 
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m Ma Mk Mr^ Fq QA D F C Ka Qm 

SU(2) chiral limit choice 880 1152 140 - 87 1.27 -- 2.85 3.16 

physical choice 938.9 1232 138.04 - 92.21 1.27 -- 2.85 3.16 

SU(3) chiral limit choice 880 1152 140 496 547 87 - 0.623 0.441 -D -3.16 

physical choice 1149 1381 140 496 547 108 - 0.8 0.47 -0.85 - 3.16 


Table I: Numerical values for the hadron masses and decay constants used in the 70 calculations. All the dimensionfull values 
are given in units of MeV. The physical choice values for SU(2) were taken as in Ref. [12] — notice the difference by a factor 2 
in the Lagrangian definitions of Ha —, whereas for the chiral limit and SU(3) we followed Ref. [2S]. The value for the coupling 
qm is calculated through Eq. m 


forward spin polarizability one needs to assume conservation of the photon energy ui = to' — incoming and outgoing 
photons have the same momenta q = q ■ In the following, the Minkowski-space vectors used are: 


=p'i^ =( to , 0 , 0 , 0 ) 
e'^=( 0 ,e) 
e*^ =( 0 ,e*) 

q^ = q'^ =7,g), 


(15) 


where q and q' are the 4-momenta of the incoming and outgoing photons, e and e* their respective polarizations, while 
p and p' are the momenta of the incoming and ongoing baryons, respectively. We work in the Weyl gauge, which 
leads to the condition p • e = 0 . 

All terms containing the expression contribute to 70 , as can be seen when comparing Eq. [^with 


f^ = ia(6xr*)-(r6*). 


(16) 


Terms like ^*0 yield a contribution of —iuja{e x e *) when projected onto the baryon states. All the other expressions 
that arise can be reduced to this simple case. 

While the full set of diagrams in the spin- 1/2 sector is gauge invariant, special care has to be taken when including 
the spin-3/2 states. The diagrams which include a minimal coupling of the photon to the A(1232) would need terms 
of higher order to fully restore gauge invariance. In order to be fully gauge invariant the Lagrangian of Eq. should 
include the covariant derivative D^Aj^ as opposed to the partial derivative only. The difference between the two 

derivatives are higher-order terms. To solve this discrepancy, we follow the solution of Ref. [ 6 ], where this problem 
has already been addressed for the case of the proton. In fact, for the neutral octet baryons the diagrams of Fig. 
are fully gauge invariant. As for the charged octet baryons, this is only the case for the diagrams with charged 
mesons. Therefore, for these baryons, the strategy is to study tw o set s of diagrams separately: on the one hand, we 
have the one-particle-irreducible diagrams of Figs. |2(b7| 2 (h)| and | 2 (i)[ which are calc ulated s umm ing over all isospin 
channels; on the other hand, the missing one-particle-reducible loop diagrams of Figs. 2 (c)| to [2(g)[ are first calculated 
only for the charged meson channels. For the other channels the isospin factor is chosen such that the ratio between 
the isospins of the one-particle-reducible and one-particle-irreducible diagrams is the same as for the charged meson 
channels. When doing so, gauge invariance is insured and the restoration procedure involves higher-order terms. 


IV Results and discussion 

The numerical results for our calculations, when including nucleons, pions and A resonances only (hadrons with 
no strangeness) are given in Table [^ where a comparison with the numerical values found by other groups is also 
given. Our calculation for the isospin- 1/2 sector completely agrees with the results of Ref. [ 12 ]. For completion, we 
also included how the 70 values vary when taking the chiral limit, where the masses were set to the best-fit chiral 
masses. We compare our results with the HBChPT results from Refs. [nun]. The discrepancy between the results 
does not lie in the parameter choice but in the heavy-mass expansion one assumes for HBChPT. 

As for the isospin-3/2 sector, our results differ from those of Ref. [ 12 ]. The reason for this is that we use a different 
counting scheme, and therefore a different set of diagrams. We also have a different Lagrangian, which directly 
sorts out the spurious spin- 1/2 contributions of the Rarita-Schwinger spin-3/2 spinor. In Ref. [ 8 ] the A(1232) was 
introduced in the same way as in the present work. One should remark that there a tree-level diagram of order 
was included, which we left out here for consistency. Without this diagram, the numerical results in Ref. [ 8 ] are 
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Figure 1: Diagrams contributing to 70 with isospin-1/2 intermediate states. The crossed diagrams are obtained by the substi¬ 
tutions uj ^ —uj and / o /*. 


in perfect agreement with ours. The decomposition of our results for the nucleon polarizabilities of Table [TT] into 
their individual parts is listed in Table m The main correction to the polarizability results comes from the tree-level 
diagrams with virtual spin-3/2 baryons, while their loop diagrams give only a small contribution. 

We extended the calculations to the SU(3) sector, again distinguishing between the case where the isospin-3/2 
resonances were included and where only octet baryons were taken into account as intermediate states. Here we 
also considered both cases: when taking the physical average values and when choosing the chiral limit, see Table [I[ 
We also obtain predictions for the hyperon forward spin polarizabilities. A full listing of the results for the octet 
baryons is given in Table [l^ and the decomposition of the results for the nucleons in Table |III[ When extending the 


model from SU( 2 ) to SU(3), one takes into account additional virtual states and different values for the parameters, 
whose impact we discuss in more detail below. Another interesting feature in SU(3) is also the appearance of the 
SU(3)-flavour forbidden photon transitions of the negatively charged octet baryons to those of the decuplet. The 
results of Table are also compared to HBChPT results (for preliminary results see Ref. [TO] and for a complete 
and improved analysis see Ref. [H]). For HBChPT the nucleon values for 70 change only slightly when going from 
SU( 2 ) to SU(3), the results remain large and positive. When changing to the covariant version (without the decuplet 
contribution) the SU(3) case leads to a reduction of the 70 results which still are positive. An additional inclusion of 
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Figure 2: Diagrams contributing to 70 with isospin-3/2 intermediate states. The crossed diagrams are obtained by the substi¬ 
tutions uj ^ —uj and / o /*. Except for the tree diagram, which has vertices of a second-order Lagrangian, all the vertices 
that appear are couplings of lowest-order Lagrangians. 


Model 

proton 

this work [12] [11] [13] [8] 

neutron 

this work [12] [11] [13] [8] 

HBChPT 

4.4 

4.4 

without A covariant chiral limit 

2.15 

3.24 

covariant physical values 

2.07 2.07 

3.06 3.06 

HBChPT 

1.7 

1.7 

with A covariant chiral limit 

-1.59(38) 

-0.59(38) 

covariant physical values 

-0.76(28) -1.74 -1.00 

0.15(28) -0.77 

experiment [3] 

— 1.01 ± 0.08(stat) ± O.lO(syst) 


dispersion relations [5] 

-1.34 

-0.38 


Table II: Numerical values for 70 obtained in the SU(2) sector in present and in other works in units of fm^. The choice 
of the numerical values for the constants for our own results can be found in Table [I] The error in our results when including 
the A(1232) resonance arises from the uncertainty in the value of the low-energy constant qm- 


the decuplet leads in both SU( 2 ) and SU(3) cases to negative 70 values closer to the empirical value for the proton of 
(—1.01 ± 0.08(stat) ± O.lO(syst)) • 10“^fm^ presented in Ref. [2- 

It is also interesting to compare the 70 results for the nucleons to those from dispersion relation studies found in 
Ref. [5] to be 7 q = —1.34 • 10 “^fm^ and 70 = —0.38 • 10 “^fm^. The inclusion of isospin-3/2 states, while already 
having an important effect in HBChPT, leads to an even better agreement with the empirical values in the case of 
fully covariant calculations, both when taking the chiral limit as well as when taking the average of the physical 
values for the constants. In fact, the difference between these two parameter sets is of higher chiral order for the 
polarizabilities. The main source of uncertainty of our results is the constant qm, whose variation leads to an error 
estimate as shown in Table [TV| We would like to stress that the results obtained here are not subject to uncertainties 
related to renormalization schemes; for the considered order there are no divergences or power-counting breaking 
terms entering into the value of 70 . 

As already discussed above, the inclusion of virtual decuplet states is crucial for an agreement of the nucleon B^PT 
polari zab ilities with phenomenological values, which is dominantly because of the tree-level diagrams, as shown in 
Table However, this is not the case for the 11“ and baryons since the photon transitions to the corresponding 
decuplet states 11*“ and are forbidden in SU(3) symmetry, and the tree-level diagrams do not appear. Hence, 
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virtual spin- 1/2 baryons 

virtual spin-3/2 baryons — tree level 

virtual spin-3/2 baryons — loops 

total 

p SU( 2 ) 

2.15 

-3.62 

- 0.13 

-1.59 

° SU(3) 

1.53 

-3.62 

- 0.05 

-2.14 

SU( 2 ) 

3.24 

-3.62 

- 0.21 

-0.59 

SU(3) 

2.28 

-3.62 

- 0.08 

-1.43 


Table III: Decomposition of the proton and neutron polarizability results in units of 10“^ fm^ into the contributions coming 
from the different sets of diagrams, when using the chiral limit for the masses and low-energy constants. The difference in 
results when using physical values or the chiral limit can be seen as systematical uncertainty. 


Model 

used values 

p n 

E" E° 

A 


^0 

without decuplet, HBChPT 

[Tulin] 

4.69 4.53 2.77 

2.54 2.44 

2.62 

0.52 

0.68 

without decuplet, covariant 

chiral limit 

1.53 2.28 0.90 

0.89 1.60 

1.09 

0.08 

0.15 


physical values 

1.68 2.33 0.93 

0.91 1.32 

1.28 

0.15 

0.25 

with decuplet, covariant 

chiral limit 

-2.14(38) -1.43(33) -2.72(33) 0.89 0.67(9) 

-1.69(28) 0.07 

-3.51(38) 


physical values 

-1.64(33) -1.03(33) -2.30(33) 0.90 0.47(8) 

-1.25(25) 0.13 

-3.02(33) 


Table IV: Numerical values for 70 obtained in our calculations, in units of 10“^ fm^ in the SU(3) sector. The choice of the 
numerical values for the constants in the covariant case can be found in Table [I] both for the chiral limit and for the physical 
average case. As for the HBChPT limit, we cite the results in Ref. m, which were later corrected in our work in Ref. M- 
The errors in the results with the decuplet arise from the uncertainty in the value of the low-energy constant qm- 


their values for the polarizabilities change only slightly by the small loop contributions with virtual decuplet baryons. 
To study the polarizabilites in BxPT, these two baryons might therefore be better suited than the proton and neutron 
since nearly all of the uncertainties coming from the inclusion of the decuplet drop out. Experimentally, it will be 
very hard to measure their polarizabilities, but it is feasable in lattice QCD. Furthermore, we want to emphasize that 
the main differences in numerical values for the proton and neutron polarizabilities in SU( 2 ) and SU(3) come from 
the choice of the parameters in Table All contributions coming from K oi rj mesons are negligible. Choosing the 
masses and constants in SU(3) as in SU( 2 ), which are equivalent parameter sets in terms of chiral counting up to the 
order will give nearly the same results. 

The addition of contributions would be the next step to further refine the calculation. 


V Summary 

We have presented an extended calculation for the spin polarizability 70 of the baryon octet. The framework we 
choose is based on manifestly Lorentz covariant baryon chiral perturbation theory, both in the SU( 2 ) as well as SU(3) 
versions. Furthermore, we explicitly include intermediate spin 3/2 states. The novel results of the present work 
concern both the SU(3) extension of fully covariant ChPT to the order and the inclusion of the spin 3/2 decuplet 
up to order Empirical results exist only for the nucleon case, where in both versions the inclusion of explicit 

decuplet states is crucial to find an agreement between phenomenology and BxPT. In particular, it is the tree-level 
diagram with virtual decuplet baryons that gives the dominant extra contribution. This also carries over to the SU(3) 
case, where all contributions from K and 77 loops turn out to be negligible. However, in SU(3) the two baryons 
and are prominent as their photon transitions to the corresponding decuplet states are forbidden in SU(3) 
symmetry. As a result, the decuplet tree-level contributions are not present and their polarizabilities in pure BxPT 
remain nearly unchanged, meaning that also most of the uncertainties connected to the decuplet inclusion drop out. 
Since experimental polarizability measurements for these baryons are unprobable, comparisons to results from lattice 
QCD would be very interesting, as polarizabilities to chiral order are leading order predictions of BxPT. The 70 
results for the hyperons, especially the ones for and can therefore serve as a benchmark for other calculations 



















in this sector. At this point it seems a necessity to extend the present calculation to the cases of the electric and 
magnetic polarizabilities aE and 13m of the nucleon and the baryon octet, especially the 11“ and where probably 
a similar situation as above occurs with respect to the inclusion of decuplet states. 
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Appendices 


A Basic notations of ChPT 


The octet matrices of pseudoscalar mesons (/>, photons Q and baryons B are given by 


J_ 7 r 0 I J_„ 




a=l 




K- 


TT+ K+ 


\ 




and 




G 0 0 ^ 


0 -I 0 


VO 0 -i 




The decuplet states are specified as 


( J_w 0 I _La 

x/ 2 ^ + ^/ 6 ^^ 


V 


E+ 

^0 




^111 ^ ^++ ^112 ^ J_A+ , = —A° , A222 = A- 

V3 V3 

A113 ^ 1 5]*+ a^2^ = T5]*o a223 = Ts*- 

\/3 \/6 \/3 

A 133 _ ^ ^*0 A 233 _ _1_ —A 333 _ q — 

" ’ " 


The Dirac tensors 7 ^^ are defined as 


7'^" = ^[7^7"] and 7^"^ = J {[7^ 7^], 7^ • 


(Al) 


(A 2 ) 


(A3) 


(A4) 


(A5) 


B Loop integrals and dimensional regularization 

The integrals for Eigs. and are taken in d dimensions and later dimensionally regularized to the normal 4- 
dimensional Minkowski space. To calculate the crossed diagrams, the simple substitutions 


ou ^ —uj and e* e 


(Bl) 
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have to be performed. To obtain the final numerical results for the forward spin polarizability, the integrands of 
the structure constants of the ^*^-terms are expanded up to the order The coefficients of the third order of the 
expansion are then used to evaluate the integrals. 

The following loop integrals are of interest in this work: 


/ 

/ 

/ 


d^z 1 

(27r)^ — A)^ 


d^z z^z^ 

{2ttY 


d^z z^z^z^z^ 

{2ttY 


{-ir\ r(n-f) 

(47r)d/2 r(n)A”-5 

(_l)n-iir (n- I - 1 ) 

(47r)d/2 r(n)A”-5-i 2 

(-l)"i r (n - f - 2) gl^''gP<^ + gPPg<'<^ + gP<Tgt'P 

(47r)d/2 r(n)A”-5-2 4 


As a result, one has to dimensionally regularize the integrals, obtaining the expressions 


Ai(A) 

A 2 (A) 


r(i-i) 

(47r)^/^A^“i 

r( 2 -f) 

(47r)^/^A^“^ 


A 3 (A) — 


A4(A) = 


r( 3 -f) 

(47r)c^/2A^“f 

r(4-f) 

(47r)d/2A4-f 


(? “ (Y) ^ - 7b + 1 + 0(e)) 

(? “ (^) + + ^^"0 

1 

167r2A 

1 

IGtt^ A^ ’ 


(B2) 


(B3) 


(B4) 


where e = 4 — d and /i is the scale parameter set to the proton mass in this work. For regularization, the minimal 
subtraction (MS) would have to be performed in the EOMS scheme, where terms proportional to 


^ + log(47r) - 7 £; + 1 (B5) 

are subtracted. It is interesting to note that in this work no diagram had to be renormalized, as at order no 
divergent or power counting breaking terms contribute to 70 . 
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